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1 .  INTRODUCTION 

The  processing  of  targets  by  the  AN/TPS-43FCV) - 1  radar  digital  target 
extractor  (DTE)  was  recently  investigated  from  operational  data  recorded  at 
the  RAAF  No  2  Control  and  Reporting  Unit  (2CRU)  in  Darwin  during  the  air 
defence  exercise  Pitch  Black  84-l(ref.l).  An  algorithm  was  employed  during  the 
analysis  to  discriminate  DTE  plots  associated  with  aircraft  targets  from  plots 
generated  by  clutter  and  internal  system  noise. 

The  algorithm  has  proven  to  be  an  effective  and  robust  tool  for  accurately 
computing  target  trajectories  from  plot  data  in  the  presence  of  a  high  clutter 
and  noise  background.  The  technique,  which  essentially  correlates  plot 
centroids  with  periodically  sampled  target  position  and  speed  estimates  from  a 
tracking  system,  offers  potential  for  the  a-posteriori  refinement  of  tracks  in 
operational  tracking  systems. 


2.  PLOT  EXTRACTION  AND  TARGET  DISCRIMINATION 

Multiple  plots  are  extracted  by  the  DTE  from  the  radar  video  during  each 
rotation,  or  scan,  of  the  radar  antenna.  Each  plot  is  a  digital  word  of  fixed 
length  which  contains  information  relevant  to  a  single  event.  The  three 
events  which  may  cause  a  plot  to  be  extracted  include:  the  detection  of  a 
target  in  the  radar  or  IFF  video;  the  detection  of  radar  jamming  from  the  jam 
strobe  video,  or  IFF  jamming  from  the  IFF  video;  and  the  traversal  of  the 
antenna  through  magnetic  north  (the  "north-mark",  or  reference  event)  Due  to 
the  presence  of  bit  flags  and  fixed  data  in  certain  fields,  jam  plots  and 
real-time  quality  control  (RTQC)  plots,  which  respectively  arise  from  jamming 
and  north-mark  events,  are  easily  discernable  from  plots  arising  from  target 
detections.  Jam  plots  are  therefore  ignored  in  the  following  discussion, 
whilst  the  important  role  played  by  RTQC  plots  in  the  timing  of  analysis 
events  is  deferred  for  the  present.  Henceforth  the  unqualified  term  "plots" 
shall  refer  to  the  category  of  plots  which  are  the  result  of  target  detections 
and  RTQC  plots  shall  always  be  explicitly  labelled.  Three  types  of  plot  are 
recognised  according  to  the  source  of  data  (ie  sensor):  primary  plots  contain 
data  only  from  the  primary  (ie  radar)  sensor;  beacon  plots  contain  data  only 
from  the  beacon  ( ie  IFF)  receiver;  and  correlated  plots  contain  data  from  both 
radar  and  IFF  which  pertain  to  the  same  ( ie  correlated)  target. 

The  most  important  data  in  each  plot  are  the  range  and  azimuth  of  the  two 
dimensional  centre  of  mass,  or  centroid,  of  the  target.  The  plot  centroid  is 
computed  by  the  DTE  from  the  pattern  of  signal  returns  in  the  radar  or  IFF 
video.  Supplementary  information  in  each  plot  may  include  IFF  code  data, 
target  height  measured  by  the  radar,  run  length  (azimuth  extent  of  the  target) 
and  DTE  internal  status  data,  depending  upon  the  type  of  plot. 

Varying  numbers  of  "false”  plots  (ie  plots  not  associated  with  known  targets) 
were  observed  in  the  output  of  the  DTE  during  the  exercise.  Most  were 
generated  by'  the  DTE  from  highly  correlated  clutter  returns  (for  example,  rain 
and  terrain)  and  internal  system  noise  where  the  correlation  distances 
approximate  the  dimensions  of  the  sliding  window  integrator  used  for  target 
extraction . 

A  number  of  different  methods  for  discriminating  target-associated  plots  from 
false  plots  was  initially  investigated.  Because  the  run  length  distribution  of 
false  plots  significantly  overlapped  that  of  genuine  plots  (of  the  same  type), 
a  statistical  discrimination  method  was  not  able  to  provide  sufficient 
reliability.  Exact  methods  were  therefore  examined.  One  such  method  relied 
upon  detecting  and  classifying  IFF  code  data  and  was  thus  roetr’-Td  to  ka^or, 
and  correlated  plots  rtt  utoful.css  was  severely  limited  because  many 
aircraft  could  not  respond  to  IFF;  they  were  either  not  equipped  with  IFF 
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transponders  or,  in  the  case  of  military  aircraft,  the  transponders  had  been 
deliberately  disabled.  The  majority  of  the  plots  recorded  were  found  to  be 
primary  plots  and  so  a  method  which  relied  instead  upon  classifying  radar 
height  data  was  examined.  The  usefulness  of  this  method  also  proved  to  be 
limited  because  the  radar  height  computer  hardware  responsible  for  supplying 
height  validation  strobes  and  height  data  to  the  DTE  appeared  to  be 
malfunctioning  during  the  exercise. 

A  method  based  upon  matching  plot  centroids  to  known,  approximate  target 
positions  was  subsequently  examined.  Its  main  attraction  is  that  it  can 
function  with  all  plot  types.  This  method  formed  the  basis  of  the  algorithm 
which  is  described  in  detail  in  Section  3  (p  3).  Although  the  reliability  of 
the  algorithm  could  not  be  mathematically  established  or  validated  because  of 
the  lack  of  suitable  data,  the  algorithm  appeared  to  work  reliably  in  practice 
when  applied  to  the  recorded  exercise  data.  With  modif icat ion ,  the  algorithm 
has  the  potential  for  application  in  operational  tracking  systems  for  the 
refinement  of  track  positions  and  the  management  of  crossing,  splitting  and 
merging  tracks. 

The  coarse  target  position  data  required  by  the  algorithm  were  obtained  from 
track  data  recorded  at  the  Myriad  tracking  system.  Figure  1  (p  20)  illustrates 
the  data  flow  paths  between  the  radar,  DTE,  Myriad  tracking  computer  and  the 
data  recording  system.  Myriad  was  the  only  source  of  track  data  available 
during  the  exercise,  but  another  (viz  independent)  source  of  data  could  have 
been  utilised  had  it  been  available.  There  was  a  significant  benefit  in  using 
Myriad  data.  Because  the  DTE  supplied  Myriad  with  plot  data  it  was  possible 
to  measure  the  distribution  of  Myriad  tracking  errors  with  respect  to  the  plot 
centroids.  The  algorithm  provided  the  means  of  identifying  the  plots  which 
corresponded  to  the  selected  targets.  Although  the  tracking  errors  measured 
by  this  approach  are  not  absolute  (viz  with  respect  to  precise  target 
positions),  the  additional  errors  introduced  by  the  DTE  in  computing  targe-, 
centroids  are  negligible.  Moreover,  the  tracking  error  of  any  target  of 
opportunity  can  be  determined  and  therefore  a  large  sample  of  errors  for 
statistical  testing  can  be  obtained.  This  method  obviates  the  need  for 
equipping  targets  with  very  sophisticated  global  position  monitor.:. g 
equipment,  which  is  only  warranted  at  the  limits  of  the  detection  range  wh'-re 
azimuthal  errors  in  the  computation  of  the  target  centroid  become  significant . 

The  algorithm  in  fact  contains  two  distinct  and  independent  procedure-., 
referred  to  as  Phases  I  and  II.  In  Phase  I,  plots  are  selected  for  further 
processing  in  Phase  II  if  their  centroids  lie  "sufficiently"  close  to 
coarse  position  of  the  target  supplied  by  the  tracking  system.  T:.-> 
pre-processing  of  plots  provided  in  this  phase  is  not  essential  to  -.he 
operation  of  the  core  algorithm  in  Phase  II  but  it  significantly  improves  its 
computational  efficiency. 

In  Phase  II,  the  distances  between  each  plot  centroid  in  one  scan  and  <•  i  .!i 
plot  centroid  in  the  following  scan  are  enumerated.  Plots  are  then  "pi:.-"i" 
if  their  distances  lie  within  a  specified  interval,  and  are  considered  to  ■ 
"feasible"  segments  of  the  trajectory  of  the  specified  target.  The  bound-  of 
the  selection  interval  are  derived  from  the  coarse  estimates  of  the  t  r  vt 
speed  supplied  by  the  tracking  system.  A  simple  heuristic  strategy 
utilised  to  further  reduce  the  number  of  feasible  trajectory  segments  an  : 
segments  are  synthesised  into  complete  trajectories,  each  of  which  corro-:  - 
to  an  uninterrupted  sequence  of  target  detections.  It  should  be  noted  : '.it 
although  the  algorithm  is  intended  to  obtain  the  trajectory  of  the  spin.::  i 
target,  the  trajectories  of  nearby  targets,  if  present,  will  also  be  obi  !:: •  ! 
However,  this  was  only  found  to  occur  whrr.  targets  were  flvi.,g  in  forma-.  . 
when  flight  paths  intersected.  in  this  situation  the  targets  were  located  ■.  „> 
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close  together  to  be  removed  by  Phase  I.  Finally,  the  extraction  time 
corresponding  to  each  plot  centroid  which  comprises  a  trajectory  is 
determined . 


3.  DESCRIPTION  OF  THE  ALGORITHM 

In  order  to  concisely  describe  the  processes  underlying  the  operation  of  the 
algorithm,  a  mathematical  treatment  is  provided  in  preference  to  a  computer 
program  listing  or  a  flowchart.  Set-theoretic  notation  was  selected  as  the 
most  appropriate  because  the  processes  involve  the  manipulation  of  objects, 
namely  plots  and  tracks.  The  mathematical  notation  is  summarized  on  page  17. 

3.1  Data  description 

The  historical  motion  of  all  of  the  targets  traversing  the  detection  volume 
of  a  surveillance  radar  over  a  finite  interval  of  time  At0  can  be 

represented  by  an  ordered  set  t  containing  the  position  and  velocity 
estimates  recorded  from  a  tracking  system  track  file.  Dependent  upon  the 
method  used  to  record  the  data,  the  set  may,  or  may  not,  contain  a  complete 
history  of  track  updates. 

Due  to  the  presence  of  false  alarms  and  reinitialised  tracks,  the  number  of 
tracks,  N,  recorded  generally  exceeds  the  number  of  different  targets 
present  during  the  tracking  interval  At0.  Reinitialisation  involves  the 

assignment  of  a  new  track  to  a  target  whose  previous  track  was  deleted  as  a 
result  of  the  target  being  undetected  by  the  sensor  for  several  scans;  the 
incorrect  placement  of  the  track  manoeuvre  gate;  or  the  selection  of  the 
wrong  plot  when  several  plots  lie  within  the  manoeuvre  gate. 

An  individual  track  history  is  a  member  of  the  partition  To  of  t,  viz, 


4>a  e  f,  for  all  a  =  1,2 . N 


whose  elements  can  be  represented  as  ordered  triples 


<Ra,Va,Ta>  e  *  for  all  i  =  1,2,..,L 

—  l  —i  i  a  a 


where  R°  and  V**  are  the  tracking  system  estimates  of  the  position  vector  r_ 
and  instantaneous  velocity  vector  3r/3t  respectively  of  the  target 
associated  with  the  i-th  update  of  track  o  at  time  t  =  T^  e  At„.  The 
indexed  track  update  times  are  chronologically  ordered  such  that 
T°  <  Ta  <  ....  <  T?  . 


A  north-mark  event  occurs  each  time  the  radar  antenna  is  aligned  to 
magnetic  north  and  is  accompanied  by  the  extraction  of  an  RTQC  plot.  The 
sequence  of  RTQC  plots  is  assumed  to  be  uninterrupted  over  the  tracking 
interval.  If  interruption  does  occur,  the  algorithm  will  not  operate 
correctly  on  a  complete  track;  however,  a  block  of  sequential  updates 
(viz  a  segment)  of  the  track  can  be  processed  provided  an  uninterrupted 
sequence  of  RTQC  plots  occurred  during  the  interval  between  the  first  and 
last  updates  of  the  segment.  A  track  may  therefore  be  processed  in  a 
piece-meal  fashion  by  two  (or  more)  applications  of  the  algorithm.  Another 
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assumption  is  that  there  is  negligible  delay  between  the  occurrence  of  a 
north-mark  event  and  the  recording  of  an  RTQC  plot.  The  set  of  RTQC  plot 
recording  times,  A,  therefore  corresponds  to  the  set  of  north-mark  times. 

There  is  also  assumed  to  be  negligible  delay  between  the  extraction  and 
recording  of  plots.  Each  plot  contains  an  estimate  of  the  range  and 
azimuth  of  the  centroid  of  a  potential  target.  The  set  P  is  defined  to 
comprise  the  ordered  pairs  <r,t>  of  plot  centroids  (viz  position  vectors  r) 
and  the  corresponding  plot  extraction  times  (tV 

The  north-mark  and  plot  centroid  data  to  be  analysed  are  restricted  to  the 

time  interval  At  =  [  T? ,  T?1  ]  (Z.  At0  .  The  time  T?  is  an  estimate  of  the 
a  L  — 

a 

time  of  extraction  of  the  plot  that  initiated  track  a.  Note  that  it  is  not 
the  actual  time  of  extraction  of  the  plot,  but  only  a  conservative  estimate 
of  the  earliest  time  of  extraction.  The  symbol  must  also  not  be  confused 
with  a  track  update  time;  it  is  strictly  the  lower  limit  of  the  time 
interval  which  determines  the  elligibility  of  plots  for  processing  by  the 
algorithm.  The  upper  limit  of  this  interval  is  the  time  of  the  last  track 
update . 

Let  A°  =  A  Pi  At  be  the  set  of  north-mark  times  over  the  duration  of 
a 

track  history  a  and  define  the  disjoint  subsets  a“  CZ  A°  whose  elements  are 

the  north-mark  times  prior  to  the  track  update  ,.ime  (where 

i  =  1,2,. . , L  ),  viz, 

a 


A*  =  {t°  s  A“:t"  <  T“} 


A  scan  is  precisely  defined  as  an  interval  of  time  between  consecutive 
north-mark  events.  The  antenna  rotation  period,  and  hence  the  duration  of  a 

scan,  is  assumed  to  be  a  constant,  T  . 

s 

The  time  at  the  start  of  the  j-th  scan  prior  to  the  i-th  track  update 
is  then  given  by 

t°  =  Max(Aa) 
i  1  1 

and 

j-1 

t°  =  Max(Aa  -  U  {t°  })  for  all  j=2 , 3 , . . . , J . , 
ij  i  in  j  ■  •  ’  t 

J  n=l 

where  J!  =  #(A°)  and  J.  =  #(A°  -  Aa  , )  for  if  1 . 

i  l  li-l 

The  notation  Max(*)  denotes  the  largest  element  of  a  finite  set  (ie  a  set 
with  a  countable  number  of  elements)  whilst  #(•)  denotes  the  number  of 
elements  in  a  set  (ie  cardinality).  The  north-mark  timing  diagram  is 
depicted  in  figure  2  (p  21). 


5 


WSRL-TN- 11/88 


Since  P  is  the  set  or  target  position  measurements  obtained  from  the  plots 
over  the  interval  At„,  the  subset  P„  CZ  P  of  measurements  over  the  duration 

of  track,  history  a  is  given  by 


P„  =  {<r,t>  t  P: t  E  At  }  . 

~  a 


3.2  Phase  I  processing 

3.2.1  Selection  of  plots 

Since  interest  lies  only  in  those  plots  that  are  likely  to  have  been 
associated  with  the  track  history  a,  spatial  filtering  is  employed  to 
eliminate  plots  that  were  not  located  within  the  vicinity  of  the  track. 

If  the  internal  processing  delays  in  the  DTE  and  the  tracking  system  are 
minimal,  the  position  (and  time  of  detection)  of  a  target  can  be 
estimated  from  the  centroid  closest  to  the  track  coordinate  (and 
corresponding  recording  time)  of  the  plot  extracted  prior  to  the  track 
update.  The  delays  in  the  DTE  and  tracking  system  vary  according  to 
their  respective  processing  loads  but  the  average  delays  are  assumed  to 
be  negligible  compared  with  the  average  delay  in  the  recording  of  the 
track  data. 

The  track  data  available  for  analysis  were  obtain  i  by  periodically 
recording  the  entire  contents  of  the  Myriad  track  file.  Since  the 
updating  of  the  Myriad  track  file  occurs  immediately  after  each  plot  is 
received  from  the  DTE,  each  record  of  track  data  captures  only  an 
instantaneous  "snapshot"  of  the  state  of  the  tracks  in  Myriad.  The 
intervals  between  these  snapshots  were  approximately  constant  and 
asynchronously  timed  with  respect  to  north-mark  events.  The  tra-.k 
update  times  are  assumed  to  correspond  to  the  snap-shot  recording  times. 
The  delay  between  any  track  update  and  its  subsequent  recording 
therefore  varies,  but  it  is  always  less  than  one  scan  period  provided 
the  interval  between  consecutive  recordings  is  greater  than  one  scan 
period.  If  the  recording  interval  is  greater  than  one,  Out  less  than 
two  scan  periods,  the  uncertainty  in  the  recording  delay  makes  it 
impossible  to  ascertain  which  scan  contained  the  plot  responsible  for  a 
given  track  update.  In  this  case,  the  required  plot  may  be  found  either 
in  the  scan  which  was  either  completed  immediately  prior  to  the  track 
update,  or  in  the  ensuing  scan. 

Assuming  that  the  track  data  recording  interval  exceeds  two  scan 
periods,  one  or  more  scans  of  plot  data  will  be  present  between  track 
updates.  The  position  of  the  target  in  the  first  scan  prior  to  the  i-th 
track  update  can  be  estimated  with  reasonable  accuracy  from  the  track 

coordinate,  R° .  The  position  of  the  target  in  any  earlier  scan  (which 

must  commence  after  the  (i-l)-th  track  update)  can  be  estimated  with 
more  accuracy  from  the  plot  data  by  a  recursive  procedure  which  will  now 
be  described.  The  procedure  enables  the  plots  which  established  the 
track  (ie  plots  extracted  immediately  prior  to  the  time  of  the  first 
track  update)  to  be  selected,  but  it  cannot  process  plots  recorded  after 
the  last  track  update.  This  is  not  a  limitation  since  any  plots 
associated  with  the  target  after  the  last  update  are  likely  to  have  been 
associated  with  another  track  and  may  thus  be  processed  accordingly. 
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Define  a  partition  {P_}  of  P0  whose  members  contain  only  the  plot  data 

in  the  j-th  ( j=l , 2 ,  .  .  .  ,  J  )  scan  completed  prior  to  the  i-th 

(i=l,2,...,L  )  track  update  for  the  target.  Figure  3  (p  21)  illustrates 
a 

the  partitioning  of  P0. 

The  sets 


P..  =  {<r,t>  £  P„:t  £  [t°  .  t°  .)  )  for  all  j  <  J. 

ij  -  i,J  +  l  ij  i 

comprise  the  plot  data  in  the  scans  that  occurred  between  consecutive 
track  updates,  and  the  sets 


pf,j. 


=  {<r,t>  £  P„:t  £  [t°_1  ^t"  j  )  }  for  i*l. 


comprise  the  data  in  those  scans  that  included  track  updates.  Because 
each  of  the  plot  centroids  in  a  given  scan  are  unique,  the  domain 

Dom(P..)  =  {r:  <r,t>£  P..}  of  each  set  P..  may  be  totally  ordered  such 
ij  -  -  ij  ij 

that  every  pair  of  elements  r,r*  £  DomiP^/)  satisfies  the  relation: 


(i)  0<r)  <  9 ( r* ) 

or 

(ii)  | r |  <  | r * |  if  0 ( r )  =  6(r*), 


where  9(*)  is  the  angular  displacement  of  a  position  vector  from 
magnetic  north  measured  in  the  sense  of  the  antenna  rotation. 

(Note  that  the  range  Rng(P.j)  =  {t:  <r,t>£  P  ^  ^  >  of  each  set  P_  can  only 

be  partially  ordered  because  two  or  more  plots  may  have  identical 
extraction  times). 


In  each  scan  the  target  is  assumed  to  be  located  somewhere  within  a 
circular  neighbourhood  centred  on  the  estimated  target  position  obtained 
by  the  recursive  procedure  previously  described  in  this  section.  Any 
plot  whose  centroid  falls  within  the  neighbourhood  is  considered  likely 
to  be  the  required  plot  and  is  consequently  selected  for  Phase  II 
processing.  The  radius  of  the  neighbourhood  is  adjusted  in  each  scan  in 
proportion  to  the  uncertainty  in  each  target  position  estimate. 


* 

The  sets  of  plot  data  P .  .  r~  P  .  .  selected 

tj  —  £J 

are  given  by 


from  the  initial 


scans 


(j  =  D 


A 

pn 


(<r,t>  e  PiJ: |r-R“|  <  p  0  > - 


The  radius  of  the  neighbourhood  in  the  initial  scans  is  chosen  to  be  a 
fixed  value  p„,  that  is  sufficiently  large  to  ensure  that  additional 

tracking  errors  introduced  by  manual  tracking  do  not  inhibit  the 
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processing  of  the  plots.  A  10  n  mile  value  was  initially  selected  when 
testing  the  algorithm.  Once  the  algorithm  was  validated  this  value  was 
reduced  by  trial  and  error. 

In  order  to  simplify  the  notation,  the  minimum  (Euclidean  I  distance 
between  an  arbitrary  point  whose  position  vector  is  r  and  the  plot 

centroids  contained  in  the  set  P..  is  defined  by 

ij 


d. . (r) 
tj  “ 


Min({ | r-r* | : 


r*  £  Dom ( P .  . ) }  ) 
tj 


where  Min(')  denotes  the  smallest  element  in  a  finite  set. 

The  target  position  estimate  in  each  initial  scan  is  given  by: 


r .  .  = 

-i  j 


R.,  if  P .  =0  (null  set),  otherwise 
- 1  1 1 


r  £  Dom(P^j)  :  |  r-R° | =d ^ ^ ( R° ) . 


The  corresponding  times  of  the  above  estimates  are 


'  i  1 


ct  A 

T.  if  P .  =0 ,  otherwise 
l  1 1 


<r  , t>£  P . , . 
-ll  ll 


In  the  other  scans  ( j=2 ,3, . . , J  . )  the  sets  of  selected  plot  data  are 
given  by 


1J 


=  { <r , t>  £ 


P.  . 
1J 


r .  .  -r 
-i,j-l  - 


V  T  Round  (It.  .  -t)/T  ]  }, 

max  s  i  ,  j  - 1  s' 


where  RoundC*)  denotes  the  nearest  integer  approximation  to  a  real 
number . 

The  radius  of  the  neighbourhood  is  calculated  to  be  the  maximum  distance 
that  the  target  could  have  travelled  in  an  integral  number  of  scan 
periods  at  the  maximum  speed  indicated  from  the  track  history,  ie 


=  Max  |v“|  . 

ic{ 1 , 2 , . .  .  ,L  } 
a 


max 
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The  target  position  estimates  are 


r.  ,  if  P . . =0 ,  otherwise 

-x.J-1  ij 


r  e  Dom(P..)  :  |r-r.  .  ,  |=d..(r.  .  ,) 
xj  '  -x.J-1  xj  -x.J-1 


and  the  corresponding  times  of  these  estimates  are 


1 J 


t.  ,  if  P..=0,  otherwise 
x.J-1  xj 


t  :  <r  .  .  ,t>£  P  .  .  . 
-xj  xj 


3.2.2  Resequencing  of  selected  plots 

Prior  to  further  processing,  the  sets  containing  the  selected  plot  data 
must  be  mapped  into  equivalent  sets  indexed  by  scan  number;  a  process 
which  essentially  changes  the  number  of  set  subscripts  from  two  to  one 

Let  k  represent  the  chronological  sequence  number,  or  index,  of  the 

scans  contained  within  the  time  interval  it  ,  so  that  k=l  denotes  the 

a 

0 

first  scan  commencing  at  time  t=t  .  Then,  k=2  denotes  the  second 

1  * J  i 

0  0 

scan  starting  at  t=t  .  for  JX>1  (or  t=t ,  .  if  J^l),  and  so  on  up 

1  (  J  1  ■  1  .  ,  J  2 

to  k=K  -<  y  J.'i-l  which  denotes  the  last  (completed)  scan  commencing 

°  l<i<L  1 

o 

at  t=tL  .2' 

o 

The  set  of  filtered  plot  data  in  scan  k  (k=l ,2 , . . . ,K  )  is  hence  given  by 


S.  =  P .  .  ,  , 
k  i  ,  j  - 1 


where  the  index  k  is  the  value  of  the  bijection 


f  :  U  [UM1.2 _ J.})-U.2,...Ko}, 

i=  1 


k=f (k , j )=1 - j+  ^  J  •  f°r  all  j=l,2,...,J.  and  for  ail 
l<mSi  m  1 


i=l ,  2 . L  . 

a 
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3.3  Phase  II  processing 

3.3.1  Computation  of  feasible  trajectory  segments 

Due  to  random  errors  in  the  processing  and  recording  of  plot  and  trie* 
data,  the  distances  travelled  by  a  target  in  each  scan  could  not  be 
accurately  determined.  It  is  assumed  that  the  distance  traversed  : r. 
each  scan  is  randomly  distributed  with  a  mean  and  standard  deviation 

given  by  VaT  and  V°T  ,  respectively,  where  V'a  and  Va  ire  the  mean  and 

U  s  o  s  u  o 

standard  deviation  of  the  tracking  system  estimates  of  the  target  speed 
computed  over  the  duration  of  the  crack  history,  viz, 


l<i<L 


and 


V  =  / 

a 


X 


|v-i'  /La-(V“) 


l<i<L 


A  trajectory  segment  is  defined  as  an  ordered  parr  of  plot  .e.ntroids 
<r,r'>,  obtained  from  two  plots  in  consecutive  scans,  and  lepresents  a 
sequence  of  two  coordinates  that  could  form  part  of  the  target 
trajectory.  A  trajectory  segment  is  termed  "feasible"  if  the  distance 

between  the  plot  centroids  lies  within  t£V°T  of  the  mean  iterance  V°T 

os  us 

traversed  by  the  target  over  a  one  scan  period,  where  £<  3  V3  is  a 

positive  constant  which  is  empirically  determined.  J  0 

The  set  of  feasible  trajectory  segments  between  the  position  estimates 
rtS ' ^=Dom(S^)  in  scan  k  and  r ' tS ' ^+1=Dom(S^+1 )  in  scan  ■k+l)  is 

therefore  given  by. 


Gk  =  {<£>-r’>E  SVS’k+l>:  HrE'|-VuOTsl-'rv“0Ts);k=l,2...,Ka-l 


An  example  which  demonstrates  the  selection  of  feasible  trajectory 
segments  is  presented  in  figure  4  (p  22). 

3,3.2  Reduction  of  the  number  of  feasible  trajectory  segments 

The  number  of  feasible  trajectory  segments  can  be  reduced  by  making  two 
general  assumptions,  namely  that: 

(1)  unique  target  positions  are  extracted  during  each  siir..  and 

(2)  the  shortest  trajectory  segments  are  the  most  lively  to  be 
associated  with  a  target. 

The  first  assumption  is  based  upon  the  extremely  low  probability  that 
two  or  more  targets  are  located  in  the  same  DTE  resolution  :>>il  during 
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the  same  scan.  Consequently,  a  feasible  trajectory  segment  in  is  only 

considered  to  form  part  of  the  same  trajectory  as  a  feasible  trajectory 
segment  in  if  it  shares  a  common  plot  centroid  in  scan  k+1.  The 

set  of  common  plot  centroids  in  each  scan  can  be  defined  by 

A,  =  RngCG,  )  f]  Dom(G.  , )  for  k  =  1,2, - K  -1 

k  k  k+1  a 

In  the  absence  of  data  beyond  scan  K^,  each  plot  centroid  in  scan  is 
considered  likely  to  terminate  a  trajectory;  thus, 


A..  -  RngfG  )  . 

k  k  ■  i 

a  a 

The  second  assumption  permits  a  trajectory  segment  to  be  selected  from 
alternative  trajectory  segments  in  G^,  on  the  basis  of  minimum  plot 

centroid  separation  distmee.  Although  it  is  more  intuitive  to  select 
those  segments  which  best  approximate  the  mean  target  traverse  distance 

V°T  ,  the  minimum  segment  length  strategy  proved  to  be  easier  to 
p  s 

implement,  particularly  in  ranking  alternative  segments  which  have  equal 
lengths,  and  yielded  satisfactory  results  with  empirical  data. 

The  strategy  operates  as  follows  and  is  illustrated  by  a  worked  example, 
in  figure  5  (p  23).  Initially,  the  trajectory  segments  are  ranked  in 
ascending  order  of  plot  centroid  separation  distance;  those  which  have 
equal  separation  distances  are  sub-ranked  in  ascending  order  of  the 
angular  offset  of  their  second  coordinate  vector;  and  those  which  have 
equal  separation  distances  and  angular  offsets  are  further  sub-ranked  in 
ascending  order  of  the  length  of  their  second  coordinate  vector.  Then 
the  first-ranked  segment,  the  one  with  the  shortest  separation  distance, 
is  selected  and  its  plot  centroids  are  noted.  Next,  the  segment  with 
the  lowest  ranking  whose  centroids  do  not  correspond  to  the  first  is 
selected.  Then  the  segment  with  the  lowest  ranking  whose  centroids  do 
not  correspond  to  the  two  previously  selected  centroids  is  selected. 
This  iterative  selection  process  continues  until  a  total  of 
i?(A^)<  Min  {#(S^'  )  ,#(S^+j  ' ) )  trajectory  segments  are  selected. 

The  process  is  formally  described  by  letting 


\  =  {<r.E,>E  Gk’  i'E  V 

denote  the  (non-empty)  subset  of  whose  range  is  identical  to  the  set 
A^.  Let  X^  be  totally  ordered  such  that  every  distinct  pair  of  elements 

<r,r'>,<r,r*>  e  X.  satisfies  the  relation: 

_  -  -  k 

(i)  | r-r'  |  <  |  r-r* | 

or  (ii)  8 (r ' )  <  0(r*)  if  |r-r'|  =  |r-r*| 

or  (iii)  | r ' |  <  |r*|  if  9(r' )  =  9(r*)  and  |r-r'|  =  |r-r*|. 

Define  (for  each  iteration  p=l ,  2  ,  .  .  .  ,  p  [  A,  ]  )  to  be 

k  K 


subsets  of 
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X,  such  that 
k 

x£p)  =  {<r,r’>E  XkP'1);  £'£  Dom(HkP’1)l>  £'"E  Rng[H£p'1}]} 

where 

x£n  =  Xk  and  h£P)=  {Min[X^p)l>. 

Note  that  in  the  above  expression  "Min"  specifies  the  element  in  the  set 
x£p)  which  satisfies  the  ordering  relation  with  respect  to  every  other 
element  in  the  set. 

A 

The  reduced  set  of  feasible  trajectory  segments  is  therefore 

given  by 


u 

p 


3.3.3  Synthesis  of  trajectories 

Trajectories  of  finite  length,  each  representing  an  unbroken  sequence  of 
target  detections,  can  be  synthesised  from  the  individual  trajectory- 

segments  in  U^G^.  The  process  is  described  below  with  the  assistance 
of  the  worked  example  in  figure  6  (p  24). 

Only  those  plot  centroids  in  the  domain  of  G^  and  not  in  the  range  of 

A 

^ ,  are  considered  capable  of  initiating  a  trajectory.  Excluded  from 

A 

consideration  are  those  centroids  in  the  domain  of  G.  which  are  either 

k 

the  intermediate  (viz  neither  the  first  nor  last)  centroids  of  the 
trajectories  initiated  prior  to  scan  k,  or  the  centroids  which  terminate 

A 

trajectories.  Initially,  each  centroid  in  the  domain  of  Gt  is 
considered  to  initiate  a  trajectory. 

Accordingly,  the  sets  of  plot  centroids  that  initiate  trajectories  in 
each  scan  k=l ,2, . . . ,K  -1  are  specified  by: 


and 


Bi  =  DomlG!) 


=  Dom(Gk+1)-Rng(Gk)  for  Ml, 
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(kl 

For  each  plot  centroid  r0'  J  (l)  t  B^,  where  4= 1,2,...  ,1^=  there 

exists  a  number  of  other  plot  centroids 


r(k)(*)  :  <r(k).(e),r(k)(e 
■“  m  ~m- 1  "~m 


>  e  G, 


k+m 


for  all  m=l,2, 


which  define  finite  sequences  in  the  domain  of  P,  viz, 

C.  (4)  =  <r^k)(4),rik)(4), _ r*k)(t)>  for  t  =  1,2, - L,  . 

K  “  _-\t  K 

0^(4)  thus  represents  the  sequence  of  plot  Centroids  in  the  4-th 
trajectory  initiated  in  scan  k. 

The  plot  centroids  r ^k ^ ( 4  ) , r (k^ (4 )  , 

respective  scans  k,  (k+1 ),...,  (k+M^)  . 

Whenever  any  interruption  occurs  in  the  extraction  of  plots  by  the  DTE, 
multiple  sequences  of  plot  centroids  are  generated  by  the  algorithm  with 
each  sequence  being  initiated  in  chronological  order  (viz,  in  subsequent 
scans).  When  more  than  one  sequence  is  produced  in  a  scan,  it  is 
possible  that  either  clutter  or  other  targets  are  present  in  the  u<-  ir 
vicinity  of  the  target.  Further  attempts  to  localise  the  target  must 
rely  upon  the  analysis  of  IFF,  radar  height  or  overall  (viz  spatiili 
behaviour  of  the  plot  sequences  obtained. 

Because  0^(4)  is  a  function  which  maps  the  set  of  integers  {0 , 1 2 } 
onto  the  domain  of  CZ  P,  the  set  of  plot  data  which  produced  the  2-th 
trajectory  segment  commencing  in  the  k-th  scan  is  given  by 


— M 


CZ.)  are  obtained  from  the 


k  l 


EkU)  =  { <r ,  t?t  Sk:  r  t  Rng[Ck(e)]} 


4.  COMPUTER  IMPLEMENTATION 

The  algorithm  was  implemented  as  a  computer  program  in  the  Statisti.il 
Analysis  System  (SAS)  version  5  high  level  general  purpose  language.  SAS  .,  >s 
chosen  because  the  plot  and  track  data  had  already  been  translated  into  'AS 
data  bases  and  many  of  the  analysis  programs  had  been  written  in  SAS.  Tin-  SAS 
features  most  extensively  used  in  the  algorithm  include  numerical  sorting, 
record-oriented  data  structures,  direct  addressing  of  data  and  the  statist-,  il 
program  library. 


5.  DISCUSSION 


5.1  Design  considerations 


Operational 
system  and 
cons idered 
significant 


conditions  and  limitations  in  the  performance  of 
the  recording  process  were  factors  that  had  to 
in  the  design  of  an  effective,  robust  algorithm, 
aspects  are  discussed  below. 


the 

bo 

Th«. 


st 
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(1)  Manual  track  initiation  was  exclusively  employed  during  the 
exercise.  The  time  taken  by  an  operator  to  manually  establish  a  track 
depends  on  his  ability  and  experience  to  recognise  the  detection  of  a 
new  target  and  enter  the  relevant  information  into  the  computer  system. 
The  delay  may  be  considerably  longer  than  the  two  scan  period  required 
for  automatic  initiation  if  the  target  visibility  is  reduced  by  clutter 
or  if  the  operator  is  overloaded  by  the  amount  of  information  to  be 
processed.  The  penalty  imposed  by  automatic  initiation  is  a  higher 
frequency  of  false  alarms.  Manual  initiation  was  employed  because  the 
false  alarm  rate  was  considered  to  be  unacceptably  high.  In  order  to 
obtain  maximal  length  trajectories,  the  algorithm  attempts  to  locate 
plots  extracted  several  scans  prior  to  the  time  of  recording  of  the 
initial  track  coordinate. 

(2)  Both  data  from  remote  sensors  and  visually  interpreted  local  radar 
video  were  manually  entered  into  the  tracking  system  at  various  times 
during  the  exercise.  When  a  track  position  is  not  manually  updated,  the 
automatic  system  performs  the  function.  The  algorithm  must  allow  for  the 
fact  that  manual  track  updates  are  generally  less  accurate  compared  with 
the  automatic  system. 

(3)  During  the  exercise  each  of  the  plots  extracted  by  the  DTE  were 
recorded  in  sequence  with  minimal  delay.  However,  only  approximately 
one  fifth  of  the  available  track  data  were  recorded.  Track  updates  were 
delayed  by  buffering  in  the  Myriad  track  file  for  up  to  one  scan  period 
and  snapshots  of  the  track  file  were  recorded  at  approximately  one 
minute  intervals.  In  addition,  the  numerical  precision  of  the 
time-stamp  used  for  marking  records  of  plot  and  track  data  was  limited 
to  one  second. 

(4)  The  numerical  precision  of  the  recorded  track  coordinates  was 
limited  to  one  nautical  mile  as  a  result  of  the  digital  truncation  of 
the  data  prior  to  recording.  However,  no  loss  of  precision  accompanied 
the  recording  of  plot  centroids. 

(5)  Some  loss  of  accuracy  occurred  in  the  positioning  of  track 
coordinates  with  respect  to  the  plot  centroids  due  to  track  smoothing. 
Track  smoothing  is  commonly  employed  in  tracking  systems  to  reduce 
random  measurement  errors.  Ideally,  track  data  for  the  algorithm  should 
have  been  tapped  from  the  input  to  the  tracking  filter. 

5.2  Assumptions 

The  assumptions  employed  in  the  algorithm  are  listed  below.  The  values  of 
adjustable  parameters  were  optimised,  where  possible,  by  trial  and  error  to 
improve  the  operation  of  the  algorithm. 

(1)  Each  scan  period  of  the  AN/TPS-43F  radar  was  assumed  to  be 

constant,  T  =  9.375  s.  The  sense  of  the  antenna  rotation  was  clockwise 
s 

when  viewed  from  above. 

(2)  The  delays  between  the  occurrence  of  a  north-mark  event  and  the 
recording  of  a  subsequent  RTQC  plot,  or  the  extraction  of  a  primary, 
beacon  or  correlated  plot  and  its  subsequent  recording,  were  assumed  to 
be  negligible. 

(3)  Due  to  delays  in  the  initiation  of  each  track,  several  plots  may  be 
extracted  prior  to  the  recording  of  the  first  track  coordinates.  The 

algorithm  therefore  commences  the  processing  of  plots  fT?-T®)  =  67  s 

prior  to  the  recording  of  a  selected  track  history  (a). 
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(4)  The  average  delay  in  the  processing  of  plots  by  the  DTE  was 
measured  to  be  less  than  0.5  s,  but  on  some  occasions,  the  delay  varied 
up  to  the  preset  maximum  value  of  5  s.  The  operation  of  the  algorithm 
was  only  affected  when  the  processing  of  plots  was  sufficiently 
prolonged  to  cause  their  extraction  in  a  subsequent  scan.  Under  these 
circumstances,  the  algorithm  treats  the  delayed  plots  as  though  they  are 
false . 


(5)  The  maximum  delay  in  the  recording  of  a  track  update  was  assumed  to 
be  one  scan  period. 

(6)  The  plot  centroid  measured  in  the  scan  which  completed  immediately 
prior  to  the  extraction  of  a  track  coordinate  was  assumed  to  be  located 
within  a  circular  region  of  radius  p0  centred  on  the  track  coordinate. 

The  radius  had  to  be  chosen  large  enough  to  allow  for  the  additional 
tracking  errors  introduced  by  the  manual  tracking  of  targets.  A  value 
of  5  n  mile  was  chosen  by  trial  and  error  using  a  large  number  of  track 
histories  because  it  was  not  known  which  targets  were  manually  tracked. 
Smaller  values  tended  to  cause  gaps  in  some  trajectories,  particularly 
in  the  vicinity  of  sharp  turns.  Larger  values  tended  to  increase  the 
number  of  short  trajectories  observed  further  away  from  the  track  due  to 
false  alarms. 


(7)  Due  to  random  errors  in  the  processing  and  recording  of  plot  and 
track  data,  the  distance  that  a  target  travelled  in  each  scan  could  not 
be  accurately  determined.  It  was  postulated  that  the  distance  traversed 
in  each  scan  by  a  target  associated  with  track  history  u  was  randomly 

distributed  with  a  mean  of  VaT  and  standard  deviation  of  V°T  ,  where  Va 

Us  as  u 

and  v“  are  the  respective  mean  and  standard  deviation  of  the  track 
speeds  over  the  duration  of  the  track  history. 


(8)  A  trajectory  segment  was  considered  to  be  "feasible"  (ie  associated 
with  the  track  history  a)  if  the  speed  estimated  from  the  separation  of 

the  plot  centroids  in  consecutive  scans  was  bounded  by  Va±£Va . 
r  u  o 

Extensive  testing  was  performed  to  select  an  appropriate  value  for  the 
constraint  parameter  £  which  would  trade  off  overall  computational  cost 
against  performance.  $  =  3  was  chosen,  but  any  value  between  one  and 
three  appeared  to  be  satisfactory. 


(9)  It  was  assumed  that  in  each  scan  the  DTE  extracted  a  plot  for  every 
target  detected  and  that  the  plot  centroids  were  unique.  This  assumption 
breaks  down  when  two  or  more  targets  are  located  within  the  same 
resolution  cell  of  the  DTE.  Under  these  circumstances  only  one  plot  is 
extracted.  The  size  of  each  resolution  cell  was  approximately  1.1°  in 
azimuth  (equal  to  the  antenna  lateral  beam-width)  and  0.25  n  mile  in 
slant  range  (equivalent  to  three  range  cell  widths).  When  the 
resolution  problem  occurs,  it  persists  for  more  than  one  scan  as  the 
targets  shift  from  one  resolution  cell  to  another.  The  persistence  may¬ 
be  relatively  short  when  target  trajectories  cross  over  (eg  when  a 
target  changes  altitude  or  azimuth  heading),  or  long  when  targets 
manoeuvre  in  tight  formation.  Under  these  circumstances,  the  algorithm 
is  constrained  to  process  the  plot  data  as  though  only  one  target  was 
present . 


(10)  A  heuristic  strategy  was  proposed  for  reducing  the  number  of 
feasible  trajectory  segments  to  be  synthesised  into  complete 
trajectories.  The  feasible  segments  are  firstly  ranked  in  ascending 
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order  of  plot  centroid  separation  distance  (or  angular  offset  and  length 
of  second  coordinate  vectors,  if  the  separation  distances  are  equal). 
The  first-ranked  segment  is  initially  selected.  The  second  and 
subsequent  segments  are  then  selected  as  follows.  The  n-th 
(n  =  2 , 3  ,  .  .  .  ,  #(  A^] )  segment  is  selected  which  has  the  lowest  ranking  and 

whose  plot  centroids  do  not  correspond  to  the  first,  or  the 
second .... ,or  the  (n-l)-th  selected  segment.  Although  this  strategy  is 
suboptimal,  it  produced  excellent  results  and  proved  simple  to  implement 
in  a  computer  program.  It  was  also  reasonably  efficient  in  respect  of 
the  amount  of  computer  processing  required.  More  complicated  strategies 
may  further  enhance  the  performance  of  the  algorithm  but  at  the  cost  of 
extra  computer  processing  time. 

5.3  The  effect  of  noise  and  interrupted  plot  sequences 

Under  ideal  conditions,  where  there  are  no  false  plots  in  the  vicinity  of 
the  selected  target  and  a  plot  is  extracted  for  the  target  in  every  scan 
over  the  duration  of  its  track  history,  the  algorithm  produces  a  continuous 
trajectory.  Rarely  were  these  conditions  actually  observed  during  the 
exercise.  Several  trajectories  were  sometimes  produced  because  targets 
sporadically  escaped  detection  by  the  radar  and  because  of  clutter. 

In  regions  of  intense  clutter  some  of  the  false  plots  not  removed  by 
spatial  filtering  in  Phase  I  were  sufficiently  correlated  from  scan  to  scan 
to  produce  feasible  trajectory  segments.  The  lengths  of  the  resultant 
trajectories  were  generally  short  compared  with  the  target  trajectories  due 
to  the  fact  that  false  plot  correlation  distances  were  short.  It  was  found 
that  most  false  plots  could  be  removed  from  the  exercise  data  by  rejecting 
the  shortest  trajectories  (viz  those  comprising  two  plot  centroids).  This 
strategy  appeared  to  be  more  effective  at  short  range  where  target 
trajectories  were  longer  due  to  more  reliable  detection. 

Plot  sequences  were  interrupted  whenever  plots  associated  with  the  target 
were  not  extracted  by  the  DTE  due  to  signal  loss.  Under  these  conditions, 
the  algorithm  continued  to  function  properly,  albeit  producing  multiple 
trajectories,  provided  the  tracking  system  did  not  lose  lock  on  the  target. 

It  is  essential  to  the  proper  operation  of  the  algorithm  that  the  sequence 
of  RTQC  plots  is  not  interrupted  over  the  tracking  interval.  RTQC  plot 
extraction  is  substantially  independent  of  target  detection  processes  when 
the  DTE  is  not  overloaded  by  false  plots.  Delays  were  occasionally  observed 
in  the  extraction  of  RTQC  plots  under  exercise  conditions  but  were  small 
enough  to  be  ignored.  The  only  breaks  in  the  RTQC  plot  sequence  observed 
during  the  exercise  occurred  when  the  radar  itself  malfunctioned, 
simultaneously  causing  the  loss  of  capability  for  tracking  and  the 
suppression  of  all  plot  data. 

5.4  Description  of  results 

Figure  7  (pp  25  to  30)  depicts  a  representative  sample  of  the  plot 
centroids  obtained  by  the  algorithm  from  the  exercise  data.  Each  plot 
history  corresponds  to  a  single  aircraft  that  was  consistently  tracked  over 
a  period  of  no  less  than  25  scans  (4  min) .  The  track  points  employed  in 
the  derivation  of  the  plot  centroids  are  shown  in  the  diagrams  to 
facilitate  their  comparison.  Note  that  the  vertical  axes  of  the  diagrams 
are  aligned  to  magnetic  north  and  the  track  coordinates  were  recorded  at 
approximately  one  minute  intervals.  The  consistent  and  accurate  estimation 
of  the  target  locations  determined  by  the  algorithm  contrast  strongly 
against  the  corresponding  track  coordinates. 
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The  sample  histories  indicate  reliable  tracking  of  single  targets  with 
tracking  errors  generally  accentuated  in  the  vicinity  of  aircraft  turns. 
History  S125  is  an  exception  because  manual  tracking  appears  to  have  been 
employed . 

An  example  of  the  inability  of  the  tracking  system  to  resolve  nearby- 
targets  and  avoid  track  seduction  is  illustrated  by  history  S 1 1 9  in 
figure  8  (p  31).  The  trajectories  produced  by  the  algorithm  indicate  that 
three  targets  were  present.  Initially  a  track  was  established  for  a  pair 
of  fighter  aircraft  flying  in  close  formation.  The  intersecting  flight  path 
of  another  aircraft  later  seduced  the  track.  Although  the  paths  of  the 
fighters  were  almost  identical  and  crossed  in  several  places,  it  was 
possible  for  the  algorithm  to  resolve  the  two  aircraft  because  they  were 
flying  nose-to-tail  (viz  their  flight  paths  were  slightly  staggered  in 
time) . 


6.  CONCLUSION 

An  algorithm  which  enables  the  sequences  of  plots  associated  with  a  selected 
target  to  be  obtained  from  historical  records  of  radar  plot  data  has  been 
described.  The  method  involves  matching  plot  centroids  to  the  approximate 
target  coordinates  derived  from  a  tracking  system.  The  algorithm  provides 
accurate  trajectories  (viz  chronological  sequences  of  plot  centroids)  and  the 
set  of  centroid,  extraction  time  pairs  which  uniquely  identify  the  plots 
associated  with  the  target. 

The  algorithm  evolved  from  a  requirement  which  arose  during  the  analysis  of 
the  Digital  Target  Extractor  of  the  RAAF  AN/TPS-43F(V) - 1  radar  system  located 
at  2CRC,  Darwin.  Because  the  operational  data  recorded  during  exercise  Pitch 
Black  84-1  sometimes  contained  a  large  number  of  false  plots,  a  reliable 
method  of  discriminating  target-associated  plots  from  false  plots  generated  by 
system  noise  and  clutter  was  required. 

The  coarse  target  position  and  speed  estimates  required  by  the  algorithm  were 
obtained  by  periodically  sampling  the  contents  of  the  track  file  of  the  Myriad 
tracking  system.  Plot  data  for  the  Myriad  system  were  supplied  by  the  DTE  and 
provided  the  opportunity  to  extend  the  analysis  to  measuring  the  distribution 
of  Myriad  tracking  errors.  The  fact  that  the  approximate  target  trajectory 
data  were  derived  from  the  same  plots  processed  by  the  algorithm  was 
coincidental.  Another  (viz  independent)  source  of  approximate  target 
trajectory  data,  had  it  been  available,  would  have  been  acceptable. 

A  set-theoretic  mathematical  description  of  the  algorithm  has  been  provided 
for  conciseness.  The  assumptions  employed  have  been  empirically  validated  by 
the  results  obtained  from  its  application  (in  the  form  of  a  SAS  computer 
program)  on  the  exercise  data.  Further  development  of  alternative  strategies 
in  Phase  II  of  the  algorithm  may  further  improve  its  performance,  but  at  the 
expense  of  additional  computer  processing  time. 

Operational,  system  and  recording  process  limitations  and  assumptions  which 
affected  the  design  of  the  algorithm  have  been  described. 

The  effectiveness  of  the  algorithm  in  discriminating  target-associated  plots 
from  false  plots  generated  by  clutter  or  system  noise  has  been  demonstrated. 
A  representative  sample  of  computed  target  trajectories  has  been  presented  to 
illustrate  the  operation  of  the  algorithm  upon  the  exercise  data. 

The  techniques  employed  in  the  algorithm  appear  to  offer  potential  for  the 
a-posteriori  refinement  of  tracks  in  operational  tracking  systems. 
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7.  MATHEMATICAL  NOTATION 


set  or  partition 
sequence 
ordered  pair 
ordered  triple 
closed  interval 
half  open  interval 
summation 
derivative 
'subset  of' 

'proper  subset  of' 

' less  than' 

'less  than  or  equal' 

' equal ' 

'not  equal' 

' equivalent ' 

'approximately  equal' 

'belongs  to'  (a  set  or  class) 

'does  not  belong  to' 

'such  that' 

'mapping'  (of  sets) 

set  union 

set  intersection 

set  product 

length  of  a  vector 

angular  displacement  of  a  vector 

nearest  integer  approximation  to  a  real  number 

domain  of  a  binary  relation 

range  of  a  binary  relation 
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MinC) 

>lax(*) 

*(*) 

0 


smallest  element  of  a  finite  set 
largest  element  of  a  finite  set 
number  of  elements  in  a  set  (cardinality) 
null  set 
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Figure  1.  General  data  flow  diagram 


Figure  3.  The  partitioning  of  the  set  P 
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This  example  assumes:  (1)  l<k<Xa'2  and 

(2)  G^  ={<a  a ' > , <b , a  1 > ,<c , b 1 > , <c , c ' > , <d , c ' > , <e ,d ' > } 

(3)  <a ' , a"> ,<a ' ,b”> , <c ' ,b"> ,<d' ,  c"> ,<e' , c"> ,<e 1 ,  d">) 

(4)  9(b)  <  0(a) 


Thus,  Rng(G^)={a' ,b' ,c' ,d' ) ,  Dom(G^+  )=t a ' ,c '  . d '  , e ' } , 
A^=(a'  ,  c'  ,  d'  }  and  4‘(A^)=3. 


Therefore,  X^^=X^=  { <a ,  a '  >  ,  <b ,  a '  >  ,  <c  ,c '  >  ,  <d ,  c '  >  ,  <e  ,d '  > )  and  {<c,c'>), 

{ <a  ,  a  '  > ,  cb ,  a  '  »  ,  ,  d '  >  }  and  {<b,a'>), 

and  X^3)=  H^3)=  {<e,d'-). 

Hence  Gk=  H^1}  y  H^2)y  H^3)=  { <  b ,  a '  >  ,<c  ,c' >  ,<e  ,d' >} 


Figure  5.  Reduction  of  the  set  of  feasible  trajectory  segments 
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This  example  assumes:  (1)  K^=6  and 


(2)  aj.bj.cj  e  S’j 

(4) 

-3 *~3 ’-3 ’-3  £ 

s 1 
b  3 

(6) 

a, , b_  t  S’  - 

“3  -3  3 

(3)  a„.b.,r  c  * '  , 

-2  -i  -t  2 

(5) 

-4  *  -4  *  ^4  ’  d4  S 

S'4 

(7j 

V^6  E  S'& 

Then,  B  ={a  •  *?]  .c^}  , 

B2=8 

.  B3=(d3), 

B,  =8 

*4 

and 

B5=(b.}. 

Setting  rg1)(l)=a!>  1 J  (  2  ^=- 1  ’  ^  f  3  )=Cj  ,  jQ3^fl)=d3.  and  £o^<1,  =  -5’ 

the  trajectories  are  therefore  given  by: 

Cl<  1  >=<ar-2,a3' VW’  Cl(2)=<^r  W-4>1  Cj(3)=^c1,c2.c3,c4>, 

C  (l)=^d.,d,>,  and  C. ( 1 )=<b_ ,b,> 

3  ~3  -4  o  -o  -o 


Figure  6.  Synthesis  of  trajectories 
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Figure  7(Contd.)- 
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